Measurement of Strange Quark Contributions to the Nucleon's Form Factors at 2^=0.230 (GeV/c)^ 
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We report on a measurement of the parity- violating asymmetry in the scattering of longitudinally polarized 
electrons on unpolarized protons at a of 0.230 (GeV/c)^ and a scattering angle of Ge = 30° — 40°. Us- 
ing a large acceptance fast PbF2 calorimeter with a solid angle of Af2 = 0.62 sr the A4 experiment is the 
first parity violation experiment to count individual scattering events. The measured asymmetry is Aphys = 
(— 5.44±0.54stat±0.26sys) x 10^^. The Standard Model expectation assuming no strangeness contributions to 
the vector form factors is Aq = (— 6.30±0.43) x 10^*. The difference is a direct measurement of the strangeness 
contribution to the vector form factors of the proton. The extracted value is + 0.225G^ = 0.039 ± 0.034 or 
Ff + O.OOFI =0.032 ±0.028. 

PACS numbers: 12.15.-y, 11.30.Er, 13.40.Gp, 13.60.Fz, 14.20.Dh 



> 



o 
o 

CD 
I 



X 



The understanding of sea-quark degrees of freedom of the 
nucleon and the interaction of the quarks with the vacuum in 
the nonperturbative low energy regime of quantum chromo- 
dynamics (QCD) is very poor even today. Since the nucleon 
has no net strangeness, any contribution of strange quarks to 
the nucleon structure observables is of particular interest be- 
cause it would be a pure sea-quark effect. For example the 
scalar strangeness content of the nucleon that gives a contri- 
bution to the mass of the nucleon has been discussed in the 
context of the E commutator which can be related to the 7I-N 
scattering amplitude 1 1|. New measurements suggest a large 
contribution of strange quarks to the mass of the nucleon on 
the order of 50%. The interpretation of the unexpected nu- 
cleon spin content results |2| suggests a sizeable contribution. 
As w —10%, of the strange quarks to the nucleon spin. 
Estimates of the strange quark contribution to the magnetic 
and electric vector form factors predict sizeable effects ac- 
cessible to experiments ||^ ||4|. Recently, two experiments, 
SAMPLE at Bates (3\ and HAPPEX at TJNAF |6], have ex- 
plored parity-violating (PV) asymmetries on the proton and 
the deuteron in two different kinematical regions. We report 
here on a new measurement at a four momentum transfer 
of 0.230 (GeV/c)^ at the Mainzer Mikrotron accelerator facil- 
ity (MAMI) |7]. The A4 experiment at MAMI is complemen- 
tary to other experiments for two reasons. Firstly its -value 
tests models predicting an enhanced strangeness contribution 
at this point |8] and secondly, for the first time counting tech- 
niques are used in a scattering experiment measuring a PV 
asymmetry. Therefore possible systematic contributions to 
the experimental asymmetries and the associated uncertainties 
are of a different nature as compared to previous experiments, 
which use analogue integrating techniques. 
Access to the strangeness nucleon vector cuiTent matrix el- 
ements is possible in the framework of the standard model 



by a measurement of the weak vector form factors ^ of 
the proton |9]. They can be expressed in terms of the known 
nucleon electro-magnetic vector form factors G^']^ and the 
unknown strangeness contribution G| ^. The interference be- 
tween weak (Z^) and electro-magnetic (y) amplitudes leads 
to a PV asymmetry AiR{ep) in the elastic scattering cross 
section for right- and left-handed electrons {<3r and Ol re- 
spectively), which is given in the framework of the Standard 
Model llOll and can be expressed as a sum of three terms, 
ALR^ep) =Av +As+Aa, with 
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Ay represents the vector coupling on the proton vertex where 
the possible strangeness contribution has been taken out and 
has been put into A, . is a term arising only if a contri- 
bution from strangeness to the electro-magnetic vector form 
factors is present and the term Aa represents the contribution 
from the axial coupling at the proton vertex due to the neutral 
current weak axial form factor G^. The quantity a represents 
(G^2^)/(47taV2). G^ is the Fermi coupling constant as de- 
rived from muon decay, a is the fine structure constant, the 
negative square of the four momentum transfer, T = 2^/(4M^^) 
with Mp the proton mass and e = [1 +2(1 +T)tan2(9e/2)]"' 
with Qe the laboratory scattering angle of the electron. The 
electro-magnetic form factors G^'^/ ^6 taken from a recent 
para metrization (version 1, page 5) by Friedrich and Walcher 
111 ill , where we assign an experimental error of 3 % to G^ and 
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G^, of 5% to G^f, and of 10% to G^. Electro-weak radiative 
corrections are included in the factors p'^^^ and ic^^ which have 
been evaluated in the MS renormalization scheme l|l2ll. We 
use a value for s| = sin^ Qw {Mz)ms of 0.23 1 1 3( 15) Q . The 
electro-weak radiative corrections to Aa have also been esti- 
mated within the MS renormalization scheme and the radiative 
corrections given in lll4(l as well as a value of Ai = —0. 1 ± 0. 1 
are included in Ga ■ Electro-magnetic internal and external ra- 
diative corrections to the asymmetry and the effect of energy 
loss due to ionization in the target have been calculated. They 
reduce the expected asymmetry in our kinematics by about 
1.3%. In order to compare to the measured asymmetry, we 
include the effects of electro-magnetic radiative corrections 
and energy loss in the target and average Aq —Av+Aa over 
the acceptance of the detector and the target length. We ob- 
tain the expected value for the asymmetry at the averaged 
without strangeness contribution to the vector form factors of 
Ao(fi2 = 0.230(GeV /cf) = -6.30 ± 0.43 ppm. 
The PV asymmetry was measured at the MAMI accelerator 
facility in Mainz |7| using the setup of the A4 experiment 
lITsll . The polarized 854.3 MeV electrons were produced us- 
ing a strained layer GaAs crystal which is illuminated with 
circularly polarized laser light 1 16|. Average beam polariza- 
tion was about 80%. The helicity of the electron beam has 
been selected every 20.08 ms by setting the high voltage of a 
fast Pockels cell according to a pattern of four helicity states, 
either (+P -P-P + P) or {~P + P + P~P). The pattern 
was selected randomly by a pseudo-random bit generator. A 
20 ms time window enabled the histogramming in all detector 
channels and an integration circuit in the beam monitoring and 
luminosity monitoring systems. The exact window length was 
locked to the power frequency of 50 Hz in the laboratory by 
a phase locked loop. For normalization, the gate length was 
measured for each helicity. Between each 20 ms measurement 
gate, there was an 80 /js time window for the high voltage at 
the Pockels cell to be changed. The intensity I = 20 /jA of the 
electron current was stabilized to better than 5/// ~ 10^^. An 
additional A,/2-plate in the optical system was used to rotate 
small remaining linear polarization components and to control 
the helicity correlated asymmetry in the electron beam current 
to the level of < 10 ppm in each five min run. 
We have used a system of microwave resonators in order to 
monitor beam current, energy, and position in two sets of 
monitors separated by a drift space of about 7.21 m in front 
of the hydrogen target. In addition, we have used a system of 
10 feed-back loops in order to stabilize current, energy 1 17|, 
position, and angle of the beam. The polarization of the elec- 
tron beam was measured using a M0ller polarimeter which 
is located on a beam line in another experimental hall. This 
M0ller polarimeter U8il has an accuracy of 2.1 %. Due to the 
fact, that we have to interpolate between the weekly M0ller 
measurements, the uncertainty in the knowledge of the beam 
polarization increased to 4%. The 10 cm high power, high 
flow liquid hydrogen target was optimized to guarantee a high 
degree of turbulence with a Reynolds-number of R > 2x 10^ 
in the target cell in order to increase the effective heat trans- 
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FIG. 1 : An energy spectram of accepted particles from the hydrogen 
target as read directly from the hardware memory of the read out 
electronics of the lead fluoride calorimeter. For the solid black curve, 
the spectrum has been corrected for the differential nonlinearity of 
the ADC. The position of the elastic scattering peak, the threshold 
for Tt'^-production and the position of the A-resonance is indicated as 
well as the lower and upper cut position for the extraction of A'f and 
Nj; as described in the text. 



fer This new technique allowed us for the first time to avoid a 
fast modulation of the beam position of the intense cw 20 /jA 
beam and it also allowed us to stabilize the beam position on 
the target cell without target density fluctuations arising from 
boiling. The total thickness of the entrance and exit aluminum 
windows is 250 /jm. The luminosity L was monitored for 
each helicity state (R, L) during the experiment using eight 
water-Cerenkov detectors (LuMo) that detect scattered parti- 
cles symmetrically around the electron beam for small scat- 
tering angles in the range of O^. = 4° — 10°, where the PV 
asymmetry is negligible. The photomultiplier tube currents 
of these luminosity detectors are integrated during the 20 ms 
measurement period by gated integrators and then digitized by 
customized 16-bit analogue-to-digital converters (ADC). The 
same method was used for all the beam parameter signals and 
these data were stored with the helicity information as a func- 
tion of time for each 20 ms time window. A correction was 
applied for the nonlinearity of the luminosity monitor photo- 
multiplier tubes. This was measured and verified separately 
by varying the beam current from 0-23 /jA several times per 
week. From the beam current helicity pair data l'^'^ and lumi- 
nosity monitor helicity pair L^'^ data we calculated the target 
density p*'^ = iJ^-^/j^-^ for the two helicity states indepen- 
dently for the analysis of the data. 

As particle detector we developed a new type of a very fast, 
homogeneous, total absorption calorimeter consisting of in- 
dividual lead fluoride (PbF2) crystals fl^ . The material is a 
pure Cerenkov radiator and has been chosen for its fast tim- 
ing characteristics and its radiation hardness |20]. This is the 
first time this material is used in a large scale calorimeter for 
a physics experiment. The crystals are dimensioned so that 
an electron deposits 96% of its total energy in an electro- 
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magnetic shower extending over a matrix of 3 x 3 crystals. 
Together with the readout electronics this allows us a mea- 
surement of the particle energy with a resolution of 3 .9 % / a/E 
and a total dead time of 20 ns. At the time of the data taking, 
511 channels of the detector and the read out electronics were 
operational. The detector modules were located in two sectors 
covering an azimuthal angle interval A(|) of 90° symmetrically 
around the beam axis. The particle rate within the acceptance 
of this solid angle is w 50 x 10'' s^'. Due to the short dead 
time, the losses due to double hits in the calorimeter are 1 % at 
20 fjA. The signals from each cluster of 9 crystals are summed 
and integrated for 20 ns in an analogue summing and trigger- 
ing circuit and digitized by a transient 8 -bit ADC. There is 
one summation, triggering, and digitization circuit per crys- 
tal. The energy, helicity, and impact information are stored 
together in a three dimensional histogram. Figure [2 shows a 
typical energy spectrum of scattered particles from the hydro- 
gen target at an electron current of 20 /jA . It was taken during 
five minutes and is a direct output of the histogramming mem- 
ory. The elastic scattering peak is clearly isolated at the high 
end of the spectrum. 

The energy resolution Oe of each group of 9 crystals was de- 
termined by analyzing such spectra. The number of elastic 
scattered electrons is determined for each detector channel by 
integrating the number of events in an interval from 1 .6 Oe 
above pion production threshold to 2.0 Oe above the elastic 
peak in each helicity histogram. The usage of those cuts en- 
sures a clean separation between elastic scattering and pion 
production or A-excitation which has an unknown PV cross 
section asymmetry. We determined the number of elastically 
scattered electrons for each helicity state (A^^ and N^) by sum- 
ming over the inner 345 detector channels which are the cen- 
ters of a full 3x3 crystal matrix. The linearity of the PbF2 de- 
tector system with respect to particle counting rates and pos- 
sible effects due to dead time were investigated by varying the 
beam current. We calculate the raw normalized detector asym- 
metry as A„w = iN^/p'^~N^/p'-)/{Nf/p'^+N^/p^). The 
possible dilution of the measured asymmetry by background 
originating from the production of which subsequently de- 
cays into two photons where one of the photons carries al- 
most the full energy of an elastic scattered electron was es- 
timated using Monte Carlo simulations to be much less than 
1 % and is neglected here. The largest background comes from 
quasi elastic scattering at the thin aluminum entrance and exit 
windows of the target cell. We have measured the aluminum 
quasielastic event rate and calculated in static approximation 
a correction factor for the aluminum of 1 .030± 0.003 giving 
a smaller value for the corrected asymmetry. 
Corrections to false asymmetries arising from helicity corre- 
lated changes of beam parameters were applied on a run by 
run base. The analysis was based on the five min runs for 
which the counted elastic events in the PbF2 detector were 
combined with the correlated beam parameter and luminos- 
ity measurements. In the analysis we applied reasonable cuts 
in order to exclude runs where the accelerator or parts of the 
PbF2 detector system were malfunctioning. The analysis is 
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FIG. 2: The extracted experimental asymmetries are shown with 
A./2-plate in or out respectively as a function of the data sample. The 
dashed line represents the value of Aq as described in the text. 



based on a total of 7.3 X 10'' histo grams corresponding to 
4.8 X 10'^ elastic scattering events. 

We extracted an experimental asymmetry from Aexp = Araw — 
a\Ai — — ajtAy — a^^hx' — flsAy' — agAZs^.. The six — 
1...6) denote the correlation coefficients between the ob- 
served false asymmetry and the electron current asymmetry 
A/ (ai), the horizontal and vertical beam position differences 
Ax, Ay (a2, 03), the horizontal and vertical beam angle differ- 
ences Axf, Ay' ((34, as), and the beam energy difference AEe 
(ag). For the analysis, the coiTelation parameters a, were ex- 
tracted by multidimensional regression analysis from the data. 
The fl, have been calculated in addition from the geometry of 
the precisely surveyed detector geometry and the a, calculated 
using the two different methods agree very well within statis- 
tics. 

The experimental asymmetry was normalized to the elec- 
tron beam polarization to extract the physics asymmetry. 
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We have taken half of our data with a sec- 



ond A-/2-plate inserted between the laser system and the GaAs 
crystal. This reverses the polarization of the electron beam 
and allows a stringent test of the understanding of systematic 
effects. The effect of the plate can be seen in Fig.|2] the ob- 
served asymmetry extracted from the different data samples 
changes sign, which is a clear sign of parity violation if, as in 
our case, the target is unpolarized. Our measured result for the 
PV physics asymmetry in the scattering cross section of polar- 
ized electrons on unpolarized protons at an average value 
of 0.230 (GeV/c)2 isAphys = (-5.44±0.54±0.26) ppm. The 
first error represents the statistical accuracy, the second er- 
ror the systematical uncertainties including beam polariza- 
tion. The absolute accuracy of the experiment represents the 
most accurate measurement of a PV asymmetry in the elastic 
scattering of longitudinally polarized electrons on unpolarized 
protons. TableUgives an overview of the applied corrections 
and the contributions to the systematical error of the experi- 
mental asymmetry. 

The interpretation of the measurement in terms of strangeness 
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TABLE I: Overview of the applied corrections and the sources of the 
experimental error in the measured asymmetry. 





correction [ppm] 


error [ppm] 


Statistics 




0.54 


Target density, luminosity 


0.58 


0.09 


Target density, beam current 


0.00 


0.04 


Nonlinearity of LuMo 


0.30 


0.04 


Dead time correction 


-0.11 


0.08 


Ai 


0.64 


0.04 




-0.05 


0.02 


Ax, Ay 
Ax', Ay' 


-0.03 


0.02 


0.03 


0.03 


Aluminum windows (H2 target) 


0.16 


0.02 


Dilution from Jt" decay 


0.00 


0.06 


Pe measurement 


-1.07 


0.11 


interpolation 


0.00 


0.19 


Systematic error 




0.26 



contribution is possible by comparing the measured physics 
asymmetry Aphys with the averaged theoretical value with- 
out strangeness contribution to the vector form factors Aq. 
The difference Aphys is proportional to an averaged com- 
bination of the Sachs form factors and the extracted value 
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0.039 ± 0.034. If one uses the Dirac 
and Pauli form factors instead, the extracted value is Fj^ + 
O.UOFj = 0.032 ± 0.028. The solid line in Fig. H illus- 
trates the possible combinations of and given by 
our result on the measured combination of + 0.225GJi^ at 
= 0.230 (GeV/c)^. The hatched area represents the eiTor 
on G;^ + 0.225G^. The measured combination is small and 
is 1 .2 standard deviations away from zero, which clearly rules 
out the pole fit type of theoretical models on the strangeness 
in the nucleon ll2li l22l . From the published result on the 
measured asymmetry of the HAPPEX collaboration |6] at 
g2=0.477 (GeV/c)2 and 9,; = 12.3°, we recalculated the com- 
bination usin g o ur parametrization for the electro-magnetic 
form factors El and yield G| +0.395G^ = 0.034 ±0.026 
and + 0.186F| = 0.024 ± 0.019. In lack of more de- 
tailed information we make the ad hoc assumption, that f| 
in the 2 -range between 0.1 (GeV/c)^ and 0.5 (GeV/cf can 
be approximated by Fj ~ —0.150 ± 0.150, coiTesponding to 
= —0.099 ± 0.154. This assumption is guided by the 
fact that this value covers within two standard deviations 
all theoretical estimates as well as the SAMPLE result from 
Isl- This yields a value for G| evaluated from our mea- 
surement of Gl-iQ^- = 0.230 (GeV/c)2) = 0.061 ±0.035 and 
FiiQ^ = 0.230 (GeV/c)2) = 0.052 ± 0.034. If one makes the 
further approximation of neglecting the -dependence in G^, 
we can combine our result with the recalculated HAPPEX re- 
sult by calculating the weighted average and yield an estimate 
of F/(0.1 (GeV/c)2 <Q^ < 0.5 (GeV/c)^) = 0.052 ±0.024 
and G|(0.1 (GeV/c)^ <Q^< 0.5 (GeV/c)^) = 0.066±0.026. 
The significance level of 2.2 a for F^ and 2.5 a for G^ which 
we reach using the assumption described above, leads us to 
the conclusion, that the combination of our measurements pre- 
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FIG. 3: 



The solid line represents all possible combinations of 

,2 



+ 0.225G^ as extracted from the work presented here at a Q' 
of 0.230 (GeV/c)^. The densely hatched region represents the uncer- 
tainty. The recalculated result from the HAPPEX published asymme- 
try at of 0.477 (GeV/c)^ is indicated by the dashed line, the less 
densely hatched area represents the associated error of the HAPPEX 
result. 



sented here with the earUer work of the HAPPEX collabora- 
tion shows for the first time evidence for the observation of a 
contribution of the strange quarks to the electric vector form 
factor of the nucleon. 
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